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Based on the Lyapunov control, we present a scheme to realize state transfer with high fidelity 
by only modulating the boundary spins in a quantum spin-1/2 chain. Recall that the conventional 
transmission protocols aim at non-stationary state (or information) transfer from the first spin to the 
end spin at a fixed time. The present scheme possesses the following advantages. First, the scheme 
does not require precise manipulations of the control time. Second, it is robust against uncertainties 
in the initial states and fluctuations in the control fields. Third, the controls are exerted only on the 
boundary sites of the chain. It works for variable spin-1/2 chains with different periodic structures 
and has well scalability. The feasibility to replace the control fields by square pules is explored, 
which simplifies the realization in experiments. 
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I. INTRODUCTION 

Quantum information processing (QIP) has been ex¬ 
tensively studied in the past decades. One of the main 
challenges in actual physical implementations of QIP is 
to transfer quantum information between different ele¬ 
ments in quantum networks. In this respect, spin chain 
systems as quantum channels [ll-Q to transfer information 
play an important role in QIP. It is believed that almost 
all spins would participate in the dynamics of an unmod¬ 
ulated spin chain Q due to spin-spin couplings, leading 
to dispersions of the chain that are harmful for quantum 
state transfer. Recently, several strategies have been pro¬ 
posed to avoid such dispersions (see, e.g., [Ml). The 
first approach is mainly based on modulating coupling 
strengths between nearest neighbors to reduce the effect 
of dispersion EMI. The second depends on the weak 
interaction of boundary spins to the remainders (bulk 
spins), making the bulk spins almost un-excited during 
the time evolution. The dynamics in the second method 
can be viewed as an effective Rabi oscillation between the 
boundary spins [isL Irzl - fll . In the other sort of schemes, 
such as modulating the Larmor frequencies on sites [ll 
or adding external potentials to the spin chain, per¬ 
fect state transfer can be obtained only for a fixed time. 
This means that a precise control on the evolution time 
is strictly required to obtain high fidelity state transfer. 
Once there exist errors in the coherent time evolution, 
the fidelity decreases sharply. 

One important progress to take the disorders into ac¬ 
count was made in Ref. [ll, where an robust state 
transmission in random unpolarized spin chains was pro¬ 
posed. This proposal is immune to some significant types 
of disorders and does not need to manipulate individual 
spins. In contrast, by manipulating the individual spin- 
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spin couplings, Zwick and co-workers investigated 

the performance for quantum state transfer and showed 
its robustness against static perturbations. 

Lyapunov-based control technique has attracted many 
attentions due to its powerful applications in manipulat¬ 
ing quantum states in quantum systems [25l - [3^ . To ap¬ 
ply the Lyapunov control, a function of the target state 
called Lyapunov function V has to be specified, where 
the control field f(t) is designed to guarantee the Lya¬ 
punov function decreasing monotonously, i.e., V < 0. 
The merit of the Lyapunov control is that the target 
is asymptotically approached, which means the manip¬ 
ulation time is not important. In this paper, we apply 
the Lyapunov-based control technique to the quantum 
state transfer from one boundary spin to another. We 
show that by manipulating the interaction strengths be¬ 
tween the boundary spins and the bulk spins in an one¬ 
dimensional spin-1/2 chain, the quantum state transfer 
can be realized. To show the robustness of this control, 
we numerically calculate the performance under the ef¬ 
fect of external perturbations. 

The paper is organized as follows. In section m we 
present a general formalism for state transfer in a spin- 
1/2 chain by Lyapunov control. In section Hill we select a 
specific eigenstate of the spin-1/2 chain for state transfer. 
The control fields are designed and the robustness against 
fluctuations is investigated in section HVl The scalability 
of the state transfer in such a spin-1/2 chain and the 
improvement on the control fields are briefly discussed in 
section IVl Finally, we give a conclusion in section IVTl 


II. GENERAL FORMALISM FOR STATE 
TRANSFER IN SPIN-1/2 CHAIN 

We start by considering an one-dimensional spin-1/2 
chain with modulated XY interactions between nearest 
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neighbors. The free Hamiltonian can be written as, 

M-l M 

= 2 E - 2 E ^rnSU^) 

m—1 m—1 

where VLm is the Larmor frequency and Dm denotes the 
coupling strength between the m-th and (m + l)-th spins. 
Sm{i = X, y, z) is the Pauli matrix and M is the length of 
spin-1/2 chain. We label the first spin and the end spin 
by 1 and M, respectively. Clearly, = 0 where 

means that the total number of spin 
up (down) is conserved. When quantum information is 
encoded on the state of spin up and spin down in the 
spin-1/2 chain, the state transfer can be equivalent to 
information transfer from the first spin to the end spin. 

By Jordan-Wigner transformation, the spin Hamilto¬ 
nian can be mapped into a spinless fermions Hamil¬ 
tonian, 

M-l M 

Hq = ^ ^ DmiCmCm+1 +c|„_|_iCm) + ^ ^ ^mCmCm, (2) 
m—1 m—1 

where Cm (c]„) represents the annihilation (creation) 
operator of spinless fermion at site m. Due to 
= 0, One Can decompose the Hilbert 
space % into subspaces H = J2j j ®9,ch of them such as 
"Hj, possessing a fixed number of spinless fermions j. We 
then choose the subspace Hq and Hi for the state trans¬ 
fer, in which it has one excitation at most. Apparently, 
the dimension of subspace Hi is equal to the number of 
spins in the chain and we denote the state with a single 
excitation at site m by |m). The Hamiltonian Hq can 
then be written in a matrix form in the basis {|m)}. 
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where the vacuum state lO) = |00---0), and |m) = 
cL|0) = |0 • ■ • OlmO • ■ • 0) [^. In the following, we study 
the problem of state transfer by using this matrix form 
of system Hamiltonian. 

To be specific, the state transfer can be obtained by 
free evolution of the Hamiltonian Hq [1,0, [l^ . It can 

be formulated as follows: The initial state of the first 
spin is a|0) -I- /3|1), while the remaining spins are in spin 
down state, 

IV'(O)) = a|0i0---00M)+/9|li0---00M) 

= a|0)+;9|l), (4) 

where |a;p-|-|/3p = 1. After a fixed evolution time t = T, 
the system would evolve to. 


Those proposals need to control the evolution time pre¬ 
cisely to get perfect state transfer. Especially, the final 
state in those proposals is not a steady state. In the fol¬ 
lowing, by Lyapunov control, we show that in the spin- 
1/2 chain one can realize steady state transfer with high 
fidelity. The control process can be illustrated as, 

1^(0)) = a|0) + /3|1) IV-(r)) = a|0) + /3|A;),(6) 

where |A/) denotes the target state and one of eigen¬ 
states of Hamiltonian Hq at the same time. Because |0) 
is the ground state of spin-1/2 chain, and we choose the 
control Hamiltonian such that the ground state remains 
unchange in the dynamics, the state transfer mainly fo¬ 
cuses on steering the initial state |1) into the target state 

We should point out that the quantum information 
encoded in the first spin is almost perfectly transferred 
to the end spin in our system, like the other protocols 
[Ml did for state transfer. The difference is that our 
resulting state is the eigenstate |A/) rather than the state 
|M) in equation ([5]). This change has many advantages 
as we mentioned in the abstract, the price we have to pay 
is the state (or information) is not localized well at the 
end spin, leading to imperfect state transfer. This can 
be improved by designing the Hamiltonain as we show 
below. 

When adding control fields fk (t) with control Hamilto¬ 
nians Hk into the system, the time evolution of spin-1/2 
chain satisfies the following Schrodinger equation {h = 1): 


K 


i?o + ^/fc(i)iffe IV')- 




(7) 


In Lyapunov control, the target state \'4 >t) is usually cho¬ 
sen as an eigenstate of the free Hamiltonian Hq, such that 
when completing the control, the target state should be 
a steady state. 


77o|^t) = A/|^/>t), (8) 

where A/ is the eigenvalue corresponding to the eigen¬ 
state I'i/'t)- To design the control fields, a Lyapunov 
function has to be chosen. By the merit of Lyapunov 
control, we choose the Lyapunov function as follows. 


H=(V^|P|V^), (9) 


where the hermitian operator P is time independent. 
Furthermore, P should commute with the Hamiltonian 
Hq, i.e., [Ho,P] = 0. With this definition, the time 
derivative of V is 

K 

v = Y,h{tmmk.pm- (10) 


The construction of the operator P is flexible, such as 


|^(T)) = a|0i0---00M)+/3|0i0---0lM) 
= a|0) + ^|M). 


P=Pf\Xf){Xf\ + Y^P^\X^){K\, (11) 


( 5 ) 
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where |Ai) {i = is the eigenstate of Hamilto¬ 

nian Hq with eigenvalue A^. To transfer a quantum state 
successfully, we should choose pi to meet the condition 
that Pf < Pi [12 • Obviously, when the system arrives at 
the target state, V reaches its minimum. Simple algebra 
shows that the control fields fk{t) = —Ak{'tlj\i[Hk,P]\'tp) 
with Hfe > 0 can assure H < 0. As a consequence, those 
control fields would transfer the quantum state from the 
spin 1 to the spin M. 

Strictly speaking, the proposal of state transfer should 
work for any state at spin 1, i.e., a and (3 in equation 
(U) could take any value. A good measure to quantify 
the performance of the state transfer then should take 
an average over all a and /3. With this consideration, we 
introduce the following averaged fidelity F{t) to quantify 
the state transfer Q, 


o K) ^ 

where r = arg|/M.i(t)| and 

set COST = 1 hereafter since t is controllable. 


( 12 ) 
1). We 


III. TARGET STATE 


Conventionally, the target state should be an eigen¬ 
state of the free Hamiltonian, such that when finishing 
the control, the target state is stationary. However, the 
perfect goal of this proposal is to transfer a state from 
one side |1) to another |M). To this extent, our pro¬ 
posal cannot be taken as a proposal for state transfer, 
nevertheless, when |A/) is almost the same as |M), our 
proposal works. So, the purpose of this section is to 
choose the parameters such that |(M|A/)| is as high as 
possible (i.e., making the occupation of spin M approach 
1 in the eigenstate |A/)). In the following, we calculate 
the eigenstates of Hamiltonian Hq in the spin-1/2 chain 
where the Larmor frequencies flm and coupling strength 
Dm in equation (|31) are periodic function of sites with pe¬ 
riod I, i.e. rtm = Alm+i and Dm = Dm+i- To be specific, 
we set the period 1 = 3. The solutions of eigenvalues and 
eigenstates of this general linear chain with open bound¬ 
ary have been given in appendix (sJl ■ We select a specific 
eigenstate of the Hamiltonian Hq so that the target state 
satisfies the condition |(M|A/)| ~ 1. By the results given 
in the appendix, the components of this specific eigen¬ 
state u can be expressed analytically as follow, 
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(13) 


where i = l,...,n a nd the correspond ing eigenvalue 
A„ = 5 [(Hi -I- H 2 ) ± From the ex¬ 

pression of eigenstate u, it is not hard to find that the 
occupations on each site (defined as |uip) are determined 


by the coupling strengths Dk{k = 1, 2, 3) and the Larmor 
frequencies Hi and H 2 . If the conditions < 1 and 

I (A 1 < 1 both satisfied simultaneously, the 

values of |u 3 i- 2 p and [rtsi-ip decrease monotonously, 
respectively. In this case, the occupation on the spin 
M is maximum or minimum, hence it is possible to let 
this eigenstate u serve as the target state. In particular, 
f{Di) = ^ is a monotonic increasing function with 
the argument Di. Thus, the value of Di should be small 
in order to obtain high fidelity state transfer, i.e., the 
boundary spins should be weakly coupled to the bulk of 
spin- 1/2 chain. 



FIG. 1: The occupation of the end spin as a function of in¬ 
teraction strength Di and D 2 . We choose the total number 
of spins M = 29. The other parameters are in units of the 
coupling strengths D 3 = 1, and fli = 1.5, fl 2 = H 3 = 0.75. 


Instead of the quadratic couplings investigated in Ref. 
M , we here consider periodic couplings for the following 
reasons. Firstly, the system is robust against small per¬ 
turbations due to the existence of gaps; Secondly, it can 
describe topological materials like topological insulators. 
In addition, systems with periodic couplings are not rare 
in practice. Figure [T] shows the occupation of end spin M 
as a function of the coupling strengths Di and D 2 . As 
expected, the occupation of end spin M becomes large 
when Di is small. For simplicity, we choose Di = 0.15, 
D 2 = D 3 = 1, Hi = 1.5, and H 2 = H 3 = 0.75 in the 
following numerical calculation. 


IV. REALIZATION OF STATE TRANSFER BY 
LYAPUNOV CONTROL 

A. Lyapunov function 

In this section, we first investigate the free evolution 
of the system starting with initial state | 1 ) without any 
external controls, that is, Hk = 0 in equation ([7]). The 
results are presented in figured! One can find that the 
state of system stays at |1) with a high fidelity. This is a 
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merit of our spin-1/2 chain: not only the final state but 
the initial state are almost stationary when the system 
does not subject to any external controls. In other words, 
it can not realize state transfer without external controls 
in our system since the system is almost still at the initial 
state. 



FIG. 2: Free evolution, namely, the time evolution without 
any controls, (a) The fidelity = (l|e“*^“*|l) versus 

time t. (b) The fidelity /M.i(t) = versus time 

t. 


Next, we control the coupling strengths between the 
boundary spins and the bulk spins, i.e., 

Hk = clck+i + k = l,M-l. (14) 

The coefficients of the hermitian operator P are chosen 
as Pi = Xi and pf = —3. The dynamics behavior of 
fidelity F is illustrated as a function of the control time 
t in figure El It approaches a steady value about 96.84% 
when the control finishes, thus we can conclude that we 
have realized state transfer with high fidelity. 


B. Robustness of state transfer 

So far, we have demonstrated that high fidelity of state 
transfer can be obtained by Lyapunov control without 
any fluctuations in the control field and uncertainties 
in the free Hamiltonian Hq. In practice, however, the 
parameters in the free Hamiltonian Hq are not easy to 
acquire precisely, and fluctuations of the control fields 
might occur in experimental manipulations, which would 
lead to errors in the state transmission. In this section, 
we focus on this issue. We begin with analyzing robust¬ 
ness against the uncertainty in the free Hamiltonian Hq, 
which can be represented by S ■ Hq, 

Hq^Hq + S-Hq. 




t 


FIG. 3: (a) Evolution of the fidelity F versus control time t. 
We have set Ak = 1 and Pf = —3 to guarantee the condition 
Pf < Pi with Pi = Xi. (b) The control field (c) The 

control field / 2 (t). 


Here S is independent random quantities manifested in 
coupling strengths 6 • Dk or Larmor frequencies S ■ flk 
{k = 1,2,3). Since the fidelity of this system will be 
impacted by the perturbations seriously when the con¬ 
trol time takes too long, we truncate the control time 
t = 1000 in the numerical calculations while the fidelity 
can approach 92%. As we only control the boundary 
coupling strengths to achieve state transfer, we empha¬ 
size on analyzing the influence of boundary parametric 
fluctuations on the fidelity at first. 



FIG. 4: The behavior of fidelity versus disorder of parameters 
related to spin 1 and spin 2. We have fixed the control time 


t = 1000. 


Figure m shows the fidelity of state transfer as a func¬ 
tion of uncertainties in the first or second spins. Those 
perturbations have slight effects on the fidelity of state 
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FIG. 5: The behavior of fidelity versus disorder of the mag¬ 
netic field fli on spin 28 and the coupling strength between 
spins 28 and 29. We have fixed the control time t = 1000. 



FIG. 6: The behavior of fidelity versus disorder of the mag¬ 
netic field 0,2 on spin 29. We have fixed the control time 
t = 1000. 


transfer since those parameters primarily affect the con¬ 
trol field /i(t), and it only plays a significant role in the 
beginning while almost vanishes at later time, as seen 
in figure [S] Figure [5] and figure |6] show the relationship 
between the fidelity and the uncertainties related to the 
28th or 29th spins. Those parameters have much im¬ 
pact on the fidelity due to the closely connection with 
the control field / 2 (t) and the occupation of end spin 
in target state. In addition, it can be found in those 
figures that the fidelity might be even higher when the 
uncertainties are very small. The reason can be under¬ 
stood as follows. In the ideal situation, the fidelity in¬ 
creases monotonously during the time evolution. Small 
uncertainties/fluctuations might not change the sign but 
change the amplitude of control fields, leading to a slight 
oscillation in the value of fidelity. As a consequence, the 
fidelity might increase due to these uncertainties. 

This can also be used to explain why a higher fidelity 
can be reached when fluctuations exist in the control 
fields /fc(t), since the fidelity sharply depends on the sign 



FIG. 7: The behavior of fidelity versus disorder of the control 
fields. We have fixed the control time t = 1000. 

rather than the amplitude of control fields in figure 0 
We find from the inset that the influence of the control 
field / 2 (t) is more sensitive than that of the control field 
fi{t) when fluctuations are small. It can be explained 
by an observation in figure |3] that there still exists small 
control field / 2 (t) when t > 1000 while the amplitude of 
control field fi{t) almost vanishes at that time. For its 
insensitive in the fluctuations of control fields, making 
the interaction time long would be in favor of obtain¬ 
ing high fidelity of state transfer. On the other hand, 
when considering the disorder of other physical parame¬ 
ters such as coupling strengths <5 • Dk or Larmor frequency 
(5 • rifc, the distribution of eigenstates of the Hamiltonian 
Hq might change. Hence the final state might not be a 
steady state and evolves even though the control fields 
disappear (or turn off). So the fidelity of state transfer 
might deteriorate with the increasing of interaction time. 
As a consequence, one should trade off the disorders of 
different parameters to truncate an appropriate control 
time in order to obtain a relatively high fidelity. 

Now, we turn to study the effect of disorder in (the 
whole spin chain) couplings on the performance of the 
state transfer. To this end, we adopt the same con¬ 
trol fields and initial state as in the case without dis¬ 
orders, but change randomly the couplings from Dj to 
Dj = {l + €)Dj, where Dj is the original couplings in the 
above analysis and j is spin site index, e is a randomly 
chosen number in the interval [-0.05, 0.05]. In other 
words, we explore the effect of disorders by re-simulating 
the state transfer with instead of Dj , and Dj changes 
randomly at each point of time for a randomly chosen 
spin site. We simulate n disorders (n = 1,2,5,20) ex¬ 
isting simultaneously in the system and show the results 
in figure [5][a) . One can discover that the scheme is ro¬ 
bust against disorders in the couplings. The disorders 
in the on-site energy would have similar effects on the 
performance. 

An interesting observations is that more disorders 
might benefit the performance of the state transfer. This 
might relate to the random number e being created in 
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No.j, J=1,...,29 


FIG. 8: (Color online) (a)Fidelity as a function of time. The 
dot line denotes the fidelity in the case without disorders, 
while the solid lines show the fidelity with disorders. Each 
line is an average over 100 simulations, where each simulation 
is obtain by introducing no disorder, 1, 2, 5 and 20 disorders. 
(b)The eigenvalues of Hamiltonian when the total spins M = 
29. The inset shows the eigenvalues of total spins M = 209. 


[-0.05, 0.05], i.e., the average of e is closely to zero. The 
physics behind the robustness of the scheme is that small 
random disorders in the couplings cannot close the gaps 
in the system since it has nontrivial topological property 
[s^, as shown in figure [8][b) . As a result, the boundary 
states appear in the system, making our control protocol 
robust, which is the essential reason why we choose the 
periodical system. 


V. SCALABILITY AND DISCUSSION 

Until now, we have just considered a chain with M = 
29 spins. A natural question arises that whether the in¬ 
teraction time increases with the increasing of the spin 
number. Since the interaction time is also closely related 
to the amplitude of control fields fk{t), we set = 1 
in the expression of fk{t) and study this problem by nu¬ 
merical calculation. Figure [HI demonstrates the relation 
between the interaction time and different total number 
of spins while the fidelity of state transfer reaches 0.93. 
It is observed that the interaction time increases approx¬ 
imately linearly with the total number of spins. 

Although the general formulas of exact diagonaliza- 
tion of matrix with period (Z > 3) in Larmor frequencies 
and coupling strengths is given in Ref. [s^, the analyti¬ 
cal expressions are complicated. For this reason, we will 
employ the numerical calculation to investigate whether 
the Lyapunov control is suited to the case for the period 
I > 3. Numerical simulations show that, it is possible 
to obtain high fidelity state transfer in a spin-1/2 chain 
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FIG. 9: The interaction time t versus different total number 
of spins M. The interaction time finishes when the fidelity of 
state transfer reaches 0.93. All the other parameters are the 
same with £gure[3l 



FIG. 10: Evolution of the fidelity F versus control time t. 
The parameters (the pink solid line) are = Us = 1.5, Q 2 = 
U 3 = U 4 = 0.75, Di = 0.15, D 3 = 0.5, D 2 = D 4 = Ds = 1. 
The parameters of the cyan dash line are the same to figure 
[ 3 ] except for the total number of spins M = 59. 


with 47 sites of period 1 = 5, and the fidelity can ap¬ 
proach 94.97% when completing the control, as shown in 
figure [ini The another feature in figure [ini is that the dy¬ 
namics behaviors of the fidelity with M = 47 and I = 5 
almost coincides with the case of M = 59 and 1 = 3, indi¬ 
cating that the interaction time does not linearly increase 
with the total number of spins in different periodic struc¬ 
ture of spin-1/2 chains. As the Lyapunov function, i.e., 
V = p/|('0|A/)p + represents the weighted 

average between the state of system and distinct eigen¬ 
states of Hamiltonian Hq. The variation of Lyapunov 
function reflects that the control fields mainly modulate 
the state transition between the eigenstates of Hamilto¬ 
nian Hq. Thus the interaction time is connected with the 
characteristic spectrum of Hamiltonian Hq rather than 
the total number of spins. Once the free Hamiltonian Hq 
of spin-1/2 chain contains an eigenstate whose occupa¬ 
tion mainly locates at the end spin, there should exist 
a Lyapunov function for the control fields to control the 
state transition between eigenstates of Hamiltonian Hq. 
As a result, it can realize state transfer by the Lyapnouv 
control, of course, the fidelity of state transmission must 
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be closely linked to the characteristic spectrum of Hamil¬ 
tonian Hq. 



t 


FIG. 11: All parameters and the expression for control func¬ 
tion are the same as in figure [3] except for the control Hamil¬ 
tonian. 


In a practical situation, the control Hamiltonians 
which exchange the interaction between the boundary 
spins and its neighbors in equation (1141) might be difficult 
to obtain. An alternative way of realizing state trans¬ 
fer is to control the Larmor frequencies at the bound¬ 
ary spins, namely, the control Hamiltonians are chosen 
as Hi = c\ci and H 2 = Figure [TT] shows the 

performance of hdelity when controlling the Larmor fre¬ 
quencies at the boundary spins. Indeed, it also can real¬ 
ize state transfer. Furthermore, we can replace the fast 
time-varying control fields by square wave pulses, which 
reduces the difficulty in experiments. The square wave 
pulses are given by 


hit) = { 


fk{t) > 0, 
fk{t) < 0. 


(15) 


By these square wave pulses, in figure [121 we show that 
the target can be reached with high fidelity as well. 


VI. CONCLUSION 

In conclusion, by Lyapunov control technique, we have 
proposed a scheme for state transfer in the spin-1/2 
chain. The controls are exerted on the couplings be¬ 
tween the boundary spins and its neighbors. It also works 
by manipulating the Larmor frequencies of the bound¬ 
ary spins. The difference between our proposal and the 
others is that the resulting state here is an eigenstate 
of the system Hamiltonian, hence it is stationary and it 
does not require to finish the transfer at a fixed time. 



t 

FIG. 12: All physical parameters are the same as in figure ITT] 
except for the control fields. We set T/ = 0.05. The control 
is terminated when the fidelity arrived at 0.96. 


In addition, the proposal is robust against many types 
of perturbations and disorders since those perturbations 
could not close the gaps of the system. This proposal 
can also be generalized to a structure with periodicity 
I > 3, provided that the chain has an eigenstate localized 
at the boundary spin (the receiver). For chains which 
have no eigenstate perfectly localized at the boundary 
spin, a high fidelity state transfer can be obtained as 
well by elaborately designing the parameters of the spin- 
1/2 chain. Finally, we have shown that replacing the fast 
time-varying control fields with square wave pulses are 
possible, which simplifies the experimental realization. 
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APPENDIX 


The diagonalization of a matrix representing the 
Hamiltonian of spin-1/2 chain with total spins M {M = 
nl -\- d), Larmor frequencies coupling strength Dm, 
periodicity I, I > d > 0 and I > 1, has been given in Ref. 























































. When I = 3 and d = 2, the matrix reads, 
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The 3n-distinct eigenvalues of this Hamiltonian can 
be obtained by solving following equation, 


(At, — Hi)(A„ — 02 )(Ai, — Hs) — (A„ — fl 2 )D^ 
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0 £>2 ••• 0 0 


\ 0 0 ■■■ D 2 0 J 


f 0 D 3 0 0 \ 

0 0 £>3 0 

\ 0 0 0 ■ • ■ 1^3 y 


1:17) 


Here [x] represents the largest integer less than or equal 
to X, the superscript “T” denotes transposition and j = 
1, 2,..., 3n. &2 and S >3 are n x (n + 1) matrices. 

The other two eigenvalues A„ satisfy the following 
equation, 

(A,-Hi)(A„-H 2 )-£)? = 0 . (18) 

The corresponding eigenvector Uj can be expressed 

Di 


U3i-ij 

U3i,j 

U3i+l,j 


A^ H 2 

0 , 

D 1 D 2 


U3i-2,j, 


{Xy — ^ 1 , 2)03 
where i = 1 ,..., n and j = 3n + 1, 3n + 2. 
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